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ABSTRACT 

We report successful fitting of a Roche model, with a surface temperature gradient following the 
von Zeipel gravity darkening law, to observations of Altair made with the Navy Prototype Optical 
Interferometer. We confirm the claim by Ohishi, Nordgren, & Hutter that Altair displays an asym- 
metric intensity distribution due to rotation, the first such detection in an isolated star. Instrumental 
effects due to the high visible flux of this first magnitude star appear to be the limiting factor in the 
accuracy of this fit, which nevertheless indicates that Altair is rotating at 0.90 ± 0.02 of its breakup 
(angular) velocity. Our results are consistent with the apparent oblateness found by van Belle et al. 
and show that the true oblateness is significantly larger owing to an inclination of the rotational axis 
of ~ 64° to the line of sight. Of particular interest, we conclude that instead of being substantially 
evolved as indicated by its classification, A7VI-V, Altair is only barely off the ZAMS and represents 
a good example of the difficulties rotation can introduce in the interpretation of this part of the HR 
diagram. 

Subject headings: stars: rotation — stars: imaging — stars: individual (aAql) — techniques: interfero- 
metric 



1. INTRODUCTION 

Altair (variously aAql, 53Aql, HR7557, HD 187642, 
of spectral type A7VI-V) is one of the brightest stars 
in the Northern sky, sharing membership in the "Sum- 
mer Triangle" with two other notable A stars. Unlike 
Vega and Deneb, Altair shows a rather diffuse spectrum 
which was early recognized to be due to a large projected 
rotational velocity variously estima ted at 242kms~ ^ 
ijUesugi fc Fukudall982), 217kms~^ (|Rover et al.l2002|) . 
and 200 km s"^ ()Abt fc Morrellflflfij^ . These estimates of 
its projected velocity, a lower limit for the true rotational 
velocity, are already a significant fraction of the breakup 
velocity, estimated near 400 km s"'^. 

Altair has become a significant object in understanding 
the atmospheres of main sequence stars at masses near 
but above that of the Sun. Specifically, Altair and a Cep 
are the two hottest stars showing Lya an d C II emis- 
sion, taken as indicators of a chromosphere l)Simon et al.l 
IT994 : Walt_gr_et__c^l- 1995) . The absence of these indica- 
tors at earlier spectral types is taken to mean that sig- 
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nificant convection disappears at this point on the up- 
per main sequence. We note that a Cep also has a high 
projected rotation velocity with estimates of 246 km s^^ 
lITTesugi fc Fuku dal1982D, 196 km s"^ lIR.over et a] .120021^ . 
and 180 km s"^ l|Abt fc Morre]llT99l listed. 

Affair's known high rotation rate has prompted at- 
tempts to measure the geometrical effects of its rota- 
tion over the years, starting with the Intensity Inter- 
ferometer (H anburv Brown et al.l I1974D . However, it 
was not until the near-IR observations with the Palo- 
ma r Tcstb cd Interfe romet er (PTI; [^olavita et al. 199!|) 
by Ivan Belle et alJ l|2001f) that a sig nificant flattening 
was detected. Comparison to classical (jvon Zeipel 1924^) 
Roche models showed the flattening was completely con- 
sistent with the observed projected rotation. 

Although this agreement between theory and observa- 
tion is nothing short of epochal, it is incomplete. Except 
near breakup, apparent oblateness, interpreted through 
Roche theory, displays the same degeneracy between 
equatorial velocity and tilt (inclination) as the apparent 
rotation velocity: one determines the quantity WegSini 
well, but not the two separately. Nor can one determine 
the sense of rotation, pro- or retrograde. Besides provid- 
ing a test of the flattening predicted by theory, oblateness 
measurements do yield the position angle of the angular 
momentum vector, the projection of that vector on the 
plane of the sky. 

In addition to flattening, Ivon Zeipell l)1924|) predicted 
that for moderate rotation stellar disks would display 
variable surface temperatures, hotter on the rotational 
axes and cooler at the equator. Speciflcally, if one de- 
fines a local effective gravity accounting for centrifugal 
acceleration, then the local effective temperature is re- 
lated to the effective gravity as T^ff oc 5e//, which is 
referred to as "gravity darkening" . With sufficient rota- 
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tion and at intermediate inclinations, gravity darkening 
predicts that stellar disks will display asymmetric inten- 
sity distributions. 

As we describe below, this prediction is of great inter- 
est in the field of optical interferomctry. Asymmetric in- 
tensity distributions produce significant imaginary com- 
ponents in the visibilities, usually represented as a non- 
trivial visibility phase. Recently developed techniques 
for recovering a closel y related quantity, " closure phase" 
l)Baldwin et all 119961: iBenson et al] 11997!) . are now be- 
ing applied to th e first round of stellar objects, (e.g., 
iWittkows ki et al.l l2?f0i;i . 

Although originally proposed as a follow-up on the 
oblateness observations, Altair was observed at the 
Navy Prototype Optical Interferometer (hereafter NPOI, 
lArmstrong et al . 1998) while the three beam combiner 
was in operation, allowing measurement of closure phase 
around one complete triangle. Examination of the data 
immediately revealed the intermediate phase angles, un- 
ambigiously signaling the presence of an asymmetric in- 
tensity profile ( Ohishi et _aL 2003) . 

Using a m odel consisting of a l imb-darkened disk and a 
bright SDot. lOhishi et ahl l|2004D demonstrated both the 
previously discovered oblateness and the necessity of in- 
cluding asymmetries in the intensity distribution. They 
argued that the probable interpretation was that of ro- 
tational flattening and gravity darkening. 

In the meantime we have become aware of some lim- 
itations in those data due to inadequate corrections for 
"deadtimes" in the avalanche photodiode detectors which 
affect the high signal levels from objects as bright as Al- 
tair. We therefore reconsider a subset of these data that 
is relatively immune to the detector proble ms, using a 
full i mplementation of von Zeipel's theory l)von Zeinell 
119241) for the model fitting, and redoing the reductions 
in a way that dramatically reduces noise in the bluest 
channels. We find that a Roche model rotating at 90% 
of the breakup angular velocity and inclined ~ 64° from 
pole-on fits the observations with high fidelity. 

We show that the parameter that sets the overall tem- 
perature scaling for the model, the effective temperature 
at the poles, Tp, is close to 8700 K for this model, and 
the polar surface gravity is correspondingly fairly high. 
This suggests that Altair is less evolved than one might 
naively expect from its spectral type and luminosity clas- 
sifications. 

In this model the equator is 1850 K cooler than the 
pole. Given that the model includes both polar bright- 
ening and a long equatorial swath of low intensity, this is 
a complex intensity distribution, and the agreement with 
the observations is a strong endorsement for the simple 
Ivon Zeipell 1)19241) theory. 

Below, we describe the new reductions, give a brief 
review of Roche theory and then present the fits. We 
note that the existence of large amounts of surface at 
near-solar temperatures suggests that the role of Altair 
(and probably a Cep) in defining the high temperature 
end of convection on the main sequence may need to 
be reconsidered. We also note the recent announcement 
that Altair is a lo w amplitude {S Set) pulsating star 
l|Buzasi et al.ll2005j) . which may give hope that astero- 
seismology will be able to put useful limits on any gra- 
dient of the angular velocity in the outer envelope. 



2. OBSERVATIONS 

Altair was observed on four nights, 25-27 May and 
1 June of 2001, with th e NPOI. These are the same ob- 
servations used bv Ohishi et al.1 i2QQ^ and lOhishi et alJ 
(2004); we refer the reader to those papers for a journal 
of observations and a description of the observing details, 
but we briefly reprise them here. We have focused here 
on the data set obtained May 25, 2001. This is by far the 
largest set of data, while the other data do not increase 
the range of hour angles observed in the first night. 

The observations used the Astrometric West (AW) , As- 
trometric East (AE), and West 7 (W7) stations, form- 
ing a triangle of interferometric baselines with lengths 
of 37.5m (AW-AE), 29.5m (W7-AW), and 64.4m (AE- 
W7) . The backend combined these three input beams to 
produce three output beams, with one baseline on each. 
The output beams were dispersed into 32 spectral chan- 
nels covering AA443 — 852 nm, although the bluest four 
channels (AA443 - 460 nm) of the W7-AW output were 
not functioning. 

The Altair observations were interleaved with obser- 
vations of a calibrator, C Aql (AOV), about 12° away 
on the sky. We initially estimated its diameter to be 
around 0.85 mas, with which it would have acted as 
a quite acceptable calibrator. However, as noted by 
lOhishi et al.1 1)2004|) , C Aql is a rapid rotator with values 
of 345kms~^ iUesugi fc FukudallT98l and 317 km s"^ 
ijRover et al.ll2002() reported. We adopted 325kms"\ 
raising the question of aspect dependent corrections to 
the squared visibilities and pha ses. This possibility was 
discussed bv lOhishi et al.l 1)2004(1 who concluded that the 
effects were not important at the level of the analysis they 
conducted, but who cautioned that the problem needed 
to be reconsidered if a detailed analysis of Altair was 
attempted with these observations. 

We have found a number of occasions in 2004 when 
C Aql was observed with a second calibrator, 7 Lyr 
(B9III). 7 Lyr is a relatively slow rotator 70kms~^, 
iRpver et all (20021) and a bit fainter {V ~ 3.24) which, 
coupled with slightly higher temperature leads to the ex- 
pectation of a symmetric, nearly unresolved calibrator 
for ^ Aql. This turned out to be correct, and to our 
surprise, we found that we were able to deduce a mean- 
ingful fit of a Roche model to ( Aql using the phase and 
visibility amplitude data. Since such results are rare, 
our report on Altair here is the first, we have decided 
to present a det ailed discussion of the case in a separate 
communication l|Peterson et al.l l2006') . 

We summarize in columns 6 and 7 of Table Q a pre- 
liminary set of the relevant Roche parameters for ^ Aql. 
Given the small angular diameter, the quoted errors pro- 
duce uncertainties in the calibration of Altair that are 
undetectable compared to other error sources. One prop- 
erty of note from these parameters is the near equator-on 
orientation of the spheroid. At most we would have ex- 
pected a few degrees uncertainty in the triple phases in- 
duced by the range of possible inclinations. However, the 
observations calibrated by 7 Lyr show that these phases 
are truly small, indistinguishable from zero at the ±1° 
le vel. The star i s app arently seen at nearly i ~ 90°. 

lOhishi et al.l l|2004|) used contemporaneous observa- 
tions of Veg cL cLS cL check star - one that should show 
a circular outline. Although we agree those observations 
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do seem to imply that Vega is circular, we are not in- 
clined to place much weight on the result. Vega is twice 
as bright as Altair in this wavelength range. Detector 
nonlinearities, already serious in the Altair observations 
as we describe in Sj^l ^.re overwhelming here. We sim- 
ply do not know how to interpret the observed squared 
visibilities. 

Finally, we note that since the NPOI records phases 
and visibilities over a range of wavelengths, there is little 
possibility of contamination by unknown companions, at 
least those within the 0'.'5 field of view of the siderostats. 
If such existed there would be clear, strongly modulated 
phases and visibilities. C Aql is a known member of a 
wide multiple system, ADS12026, with companions no 
closer than 5 arcseconds and magnitude differences no 
less than 8. 

3. DATA REDUCTION 
3.1. Incoherent Integration 

The NPOI observes interference fringes by modulating 
the optical path on the delay line for each array element, 
using a triangle-wave pattern at a frequency of 500 Hz. 
The resulting modulation of the intensity is detected in 
8 bins evenly spaced over one fringe in each channel by 
Avalanche Photo Diodes (APD). The phase of the inten- 
sity modulation changes on time scales of milliseconds 
since the fringe tracker employed by NPOI tracks the 
envelope of the (bandwidth limited) fringe packet rather 
than the fringe phase. 

The data in the delay bins were processed to produce 
the complex visibility and squared visibility modulus 
for each baseline at each wavelength. From these the 
triple product V123 exp i(j>ci can be calculated, where the 
triple amplitude 1^23 = |Vi||V2||Vi3| is the product of the 
amplitudes of the complex visibilities of the individual 
baselines and the closure phase (p^i — '/'i + 02 + 03 is 
the sum of the individual phases. Although the base- 
line phases themselves are affected by atmospheric tur- 
bulence, those effects cancel in the sum of three phases 
around a closed triangle, so the closure phase preserves 
information about the source structure. These data 
products are produced for each 2 ms cycle of delay mod- 
ulation. I n the standar d incoherent integration as de- 
scribed bv lHummel et a l. (1998), the squared visibilities 
and complex triple products are summed to provide av- 
erage values in one second intervals. 

3.2. Coherent Integration 

We employed a new algorithm for the coherent inte- 
gration of the complex visibilities of the NPOI first pre- 
sented by Hummel ct al. ( 2003). Compared to the in- 
coherent integration of the squared visiblities, coherent 
integration achieves a higher SNR of the averages due to 
the larger number of photons detected in a coherent sam- 
ple of the fringe. We have exploited this fact to recover 
meaningful results from all NPOI spectrometer channels, 
while the channels on the blue side of about 560 nm had 
usually been discarded in incoherent reductions due to 
the insufficient number of photons detected during a 2 ms 
instrumental integration time. For the coherent integra- 
tion time we selected 200 ms, and the resulting complex 
visibilities were both combined to form complex triple 
products and transformed individually into squared am- 
plitudes of the modulus. Every ten samples of these 



quantities were then averaged (averaging real and imag- 
inary part of the complex triple products separately) for 
a total integration time per data point of 2 s. 

The alignment of the raw visibility phasors neccessary 
before integration in order to avoid detrimental coher- 
ence losses was performed as follows. Two steps are nec- 
essary to rotate the phasors onto a common fringe in 
order to enable a phase tracking algorithm. 

First, average power spectra of the channeled visibility 
as a function of delay were computed for 10 ms intervals. 
Their maxima, corresponding to the group delays of the 
fringe packets, are not zero but have a typical RMS on 
the order of one micron as the NPOI group delay fringe 
tracker tries to center the fringe, but does not lock onto 
its phase. 

Second, from the deviation in position of the fringe 
from the estimated geometrical value, which is non-zero 
and has a typical RMS on the order of 10 /im due to 
atmospheric refractive index fluctuations, we estimated 
the differential amount of air and thus the phase shift be- 
tween the peak of the envelope of the fringe packet and 
the nearest fringe peak. In other words, this phase is the 
phase of the complex Fourier transform of the visibility 
as a function of wavenumber. The modulus of this trans- 
form peaks at the value of the group delay. The phase 
of the transform at this delay is called the group delay 
phase. 

We converted the group delay phase to a delay using 
the mean wavelength of the white-light fringe, and added 
it to the group delay. Rotation of the visibility phasors 
of different channels by an angle corresponding to the 
ratio of this delay value and the wavelength of the spec- 
trometer channel will align them on the same fringe. At 
this point, the algorithm implements a photon- noise lim- 
ited off-line fringe phase tracker enabling the use of much 
longer coherent integration times. 

3.3. Baseline Bootstrapping 

We used an important modification of the above 
procedure by applying the baseline bootstrapping 
method, a design feature of the NPOI interferometer 
ijArmstrong et al] IT9 98) . It exploits the fact that the 
sum of the fringe delays along a closed loop of baselines 
is zero (if the same fringe is identified on each baseline) . 
Therefore, if a long baseline in a multi-telescope array 
sees a low contrast fringe due to, e.g., object extension, 
and this baseline involves two telescopes which are at 
the same time involved with other telescopes of the ar- 
ray on much shorter baselines seeing much higher fringe 
contrast, the fringe delay of the long baseline can be 
computed from the fringe delays on the shorter baselines 
which are "bootstrapping" the long one. In the simple 
case of the observations described on Altair, the fringe 
delay on the long 64 m W7-AE baseline is just the dif- 
ference between the fringe delays on the shorter AE-AW 
and W7-AW baselines. 

3.4. Averaging and Editing 

The 2 s data points are edited for outliers as described 
by Hmnmel ct al. (1998). The final averaging is done 
over the full length of a pointing (called a scan at NPOI) 
which lasts typically 90s. The computatio n of the for- 
mal errors also follows lHummel et al.l l)1998j) . except that 
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we have implemented a different approach for the com- 
plex triple products based on a suggestion by D. Buscher 
(priv. comm.). Under simple assumptions, the error of 
a complex triple product is described by an error ellipse 
which has one axis aligned with the triple product pha- 
sor. Assuming this, we compute the error of the triple 
product as the error of the imaginary and real parts of 
the mean after applying a rotation of all phasors by the 
mean triple product phasor. The errors of amplitude and 
phase of the triple product are then equal to the error of 
the real part and the error of the imaginary part divided 
by the respective amplitudes. 

3.5. Detector Non-linearity 

A source of systematic error comes from deadtime in 
the pulse counting electronics controlling the avalanche 
photodiode detectors. These systems saturate at about 
1 MHz and display significant non-linearity in the appar- 
ent count rates as they near this limit. The nominal de- 
sign of the detector systems included a r = 200 ns dead- 
time, but we have subsequently found that not only do 
those time constants vary significantly, channel to chan- 
nel, but also that in a given channel they depend on the 
mean signal level because of the effects of heating. We 
believe that it will be possible to model and remove these 
effects, but some effort is involved, which we will report 
on in the future. Unfortunately, these problems, which 
do not affect the fainter objects usually observed by the 
instrument, were not recognized at the time the Altair 
observations were made. 

However, we believe that through a rather unique set 
of circumstances the phase and some of the amplitude 
measurements acquired during the 2001 observations are 
to first order free of the effects of these non-linearities. 
One reason was that during these observations only three 
stations were in use and the three spectrographs recorded 
single baseline data. 

The other reason was that the amplitude and phase 
measurements from each channel were accomplished with 
a simple Discrete Fourier Transform. By dithering the 
optical delay at a frequency of w, the signal was modu- 
lated according to 

I{t) = lo [I + V cos {ojt + (I))] (1) 

where V is the (instrumental) amplitude of the visibil- 
ity and the instantaneous phase. The detector system 
responded to the modulated signal according to 

Where Q is the quantuum efficiency, N(t) is the apparent 
pulse rate and No{t) = QI{t) is the true photon detec- 
tion rate. In the linearized form we assumed tNo << 1, 
which for a 300 kHz count rate, typical of the wider red 
channels on the higher visibility baselines, is adequate 
to 1% or better. Substituting equation ^ into 13 clear- 
ing the quadratic cosine using the half angle formula and 
collecting terms, this becomes 

N{t)^N^ [1 - riVo (1 -f UV2)] 

+NfiV{l - 2tNq) cos{ujt + (/)) 

^ cos(2wt -I- 20) (3) 



where Nq = QIq. A DFT at the dither frequency now 
extracts an amplitude different than the nominal NqV 
(and after division by the nominal mean signal produces 
an estimate that can differ significantly from the true 
visibility amplitude). However, the phase comes through 
the process unaffected, as do the frequencies of the min- 
ima in the visibility amplitudes. By ffoating the overall 
amplitudes in the data reduction (see below) we retain 
the important spatial scale information contained in the 
minima. But the most important conclusion is that the 
phases may be assumed to be essentially free of detector 
induced biases. 

3.6. Visibility Calibration 

The degradation of the measured visibilities due to at- 
mospheric and instrumental effects is measured, as with 
all interferometers, by observing calibrator stars with di- 
ameters as small as possible to reduce unc ertainties in the 
visibil ity estimates for them. As shown bv lHummel et all 
the NPOI visibility amplitudes sometimes show a 
negative correlation with the RMS of the delay line mo- 
tion which is related to the seeing. But at other times, 
instrumental effects which correlate with time or other 
systematic effects which correlate with position on the 
sky (e.g., hour angle) can dominate the visibility vari- 
ations. Therefore, formal photon-noise based visibility 
errors usually require the addition (in quadrature) of a 
calibration error which is derived from the residual vis- 
ibility variations of the calibrator after calibration. For 
the amplitude calibrations, we smoothed the calibrator 
visibilities with a 20min Gaussian kernel in hour angle, 
and obtained calibration errors ranging from about 4% at 
the red end to 15% at the blue end of the spectrometers. 
For the closure phase, we used the same smoothing tech- 
nique but applied to the calibrator phases as a function 
of time. 

During model fitting, however, this standard procedure 
yielded incon sistent fits. One can already see this from 
the results of lOhishi et~al] lj2004() (their Fig. 6), where 
amplitudes can be systematically high or low with the 
important characteristic that the deviation is very con- 
sistently independent of wavelength. The effect is exac- 
erbated by the fact that due to the brightness of Altair, 
the formal amplitude errors are quite small. The reason 
for the scan-to-scan variations is most likely the same as 
for the residual variations of the calibrator after calibra- 
tion, except that there is no perfect correlation due to 
target and calibrator not being at the same location in 
the sky. (Past experience has shown that visibilities do 
correlate quite well if the calibrator is very near the tar- 
get.) Therefore, we allowed the calibration for each scan 
and baseline to float by applying "achromatic" calibra- 
tion factors to improve the fit between data and model. 
We will discuss the implications for the model fitting in 
Section^] Finally we note that an error in the calibrator 
star diameter produces chromatic errors across the NPOI 
spectrometer as do uncompensated deadtime corrections, 
neither of which can be removed by the calibration fac- 
tors. 

4. MODELING 

4.1. Roche Spheroids 

The theory for the equilibrium shapes and surface 
properties of rotating stars was first presented 80 years 
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ago l)von ZeiDellll924|) assuming solid body rotation and 
a point source gravitational potential. This model has 
proved quite successful in describing the figures of stars 
in close binary systems, where tidal effect s to first orde r 
produce the same distortions as rotation (Collins 198^1). 

We use that model here, but note its limitations. First, 
there is no a priori reason that stars should rotate as 
solid bodies, and the surface layers of the Sun have long 
been known to rotate differentially. However, among the 
early results of helioseismology was the discovery that the 
transition from the outer convection zone to the inner 
radiative layers coin cided with an abrupt transition to 
solid body rotation ijSpiegel fc Zahnlll992D . Since early- 
type stars have radiative envelopes and relatively small 
(in radius) convective cores, one might expect solid body 
rotation to be a good approximation for the external lay- 
ers of early type stars. 

There has been some observational support for this 
expectation. IReiners fc R over (2004) have analyzed the 
rotational profiles of a large number of A stars looking 
for evidence of differential rotation following a solar-type 
latitudinal dependence. In the 78 stars for which the de- 
termination could be made, they found 4 objects where 
pecularities were seen which might be from differential 
rotation (or other causes). However, 95% of the line 
shapes were fully consistent with solid body rotation. 

In addition there is the long known consistency 
between the largest rotational velocities measured in 
the early-type stars and the predicted maximum ro- 
tation velocities associated with "equatorial breakup" 
dFrcmat ct al. 2005). In recent decades (e.g., Tassoul 
Il978ft it has been demonstrated that rotation laws other 
than rigid rotation do not generally impose maximum 
rotation velocities. 

The second limitation, one which we will spend some 
time on, involves the exponent in the Teff — geff re- 
lation (geff, the effective gr avitv, includes ce ntrifugal 
terms) . In the original work, Ivon Zeipell 1)1924(1 consid- 
ered the case of a fully radiative envelope, deriving the 
well known "gravity darkening" relation T^ff oc g^ff- 
ILucvl (|l£^2l) reconsidered the problem in fully convec- 
tive stars, deriving a much reduced gravity dependence, 
Tf,ff (X g^ ff- Othe r approximations lead to other expo- 
nents fsee IReinersll200d for refer ences). In our nom inal 
calculations we adopt the original Ivon Zeipell llT92^ pre- 
scription. As we shall show, the Altair observations bear 
significantly on this issue. 

Eve n in the limit of rigid rotation, the Ivon Zeioeil 
()1924ll theory is only first order in rotation rate. Distor- 
tions of the interior figure, allowing some gravitational 
quadrapole contribution, can be expected as rotation 
rates approach breakup. That is even more likely if there 
are significant deviations from solid body rotati on, even 
if con fined to the inner convective regions (e.g.. iTassoull 
IT978I) . 

Probably even more relevant are the effects of radia- 
tion pressure, which is treated simplistically in the the- 
ory, and stellar winds. Significant envelope extension due 
to radiation can be expected in the low effective surface 
gravity regions of rapidly rotating stars. And we would 
certainly expect a dramatic increase in mass loss at the 
equator, both due to enhanced convection as gas temper- 
atures decrease into the solar and sub-solar regime. We 



ignore all these effects here, except to acknowledge the 
limitations inherent in this theory. 

4.2. Roche Models 

We have constructed a suite of programs to evalu- 
ate the run of specific intensity across the surface of a 
Roche s pheroid. The definitions of the various angles 
are from ICoUinsI (^63). Otherwise, we follow the pre- 
scr iption for th e surface figure and the notation given 
bv iHardorp fc Strittmatter (19 6j) with one exception: 
following the discussion bv IHardorp fc Schold ^^1971^') we 
take the polar radius, as a fixed parameter. Specifi- 
cally we do not allow it to be a function of the fractional 
rotation. 

The modeling requires that we specify six quantities: 
the ratio of the angular velocity to that of breakup, 
uj ~ Vl/VlB^ the inclination (or tilt) of the rotational 
axis, i, defined such that i = is pole-on, the posi- 
tion angle, PA, of the pole on the sky (measured North 
through East), the angular diameter of the polar axis, 
9p, the effective temperature at the pole, Tp, and the 
surface gravity, or more commonly the logarithm of the 
surface gravity (cgs), at the pole, log gp. From the re- 
lations in the cited references it is then possible to cal- 
culate the radius of the star for a given stellar latitude, 
9, and hence the surface gravity geff and effective tem- 
perature {T^{9) = T^{ geff/gp)) at that latitude (see 
IHardorp fc Strittmatteilir968l for the definition of geff)- 

Note that for stars with accurate parallaxes like Altair, 
specifying the polar surface gravity and angular diame- 
ter is equivalent to specifying the mass and linear (polar) 
radius. These in turn fix the breakup angular velocity, 
ilg ~ {8/27)GM/R^, along with the equatorial and pro- 
jected velocities (when the inclination is specified). Fi- 
nally, it is useful to recall the relation between the polar 
radius and equa torial radius at breakup angula r velocity: 
R^ B = 2,Rp/2 (IHardorp fc StrittmatteJll968j) . Accord- 
ing to this first order theory the maximum rotational 
flattening is 2/3. 

The model definition is completed by specifying, for 
each wavelength and surface point, the specific intensity 
at the angle of the line of sight from the local normal. 
As noted by the early authors, plane-parallel model at- 
mospheres are entirely adequate in the context of Roche 
models for stars on and near the main sequence. The only 
exception is that these models develop a cusp at the equa- 
tor at critical rotation velocity. However, the cusp does 
not appear until fractional rotation velocities of a; = 0.99 
or larger and, as we have indicated, the breakdown of the 
plane-parallel approximation is only one of several prob- 
lems with the model in this limit. A number of auxiliary 
quantities are also calculated, both integrated over the 
spheroid and as seen from that particular direction, in 
addition to the predicted complex visibilities. An exam- 
ple of the forme r is the i ntegra ted luminosity, calculated 
as described in iCoUiii^ l)1963f) . Examples of the latter 
include the magnitudes in various bandpasses. 

In practice we need to solve the inverse problem: given 
a point on the sky, (a, 5), determine whether the point 
is on the stellar disk and if so what the corresponding 
latitude and radius are. We have solved this problem 
explicitly using simple iteration, and first and second or- 
der versions of Newton-Raphson iteration. The routines, 
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written in C, are quite flexible, reasonably fast and freely 
available (from the first author). 

The properties of Roche models for isolated rotat- 
ing early-type stars have recen tly been reviewed by 
iDomiciano de Soiiza et alT (2002), who use a slightly 
different but completely equivalent parameterization. 
Those authors summarize some of the results from their 
models for massive stars, which have provided a useful 
check on our own routines. 

4.3. Implementation of Roche code 

The Roche code consists of a library of functions writ- 
ten in C, with a main function enabling its use as stan- 
dalone software and a wrapper enabling it to be called 
from within the NPOI standard data reduction software 
OYSTER. The Roche spheroid parameters are part of 
the standard hierarchical model format of OYSTER and 
are passed to the Roche code along with pointers to 
the extensive tables of linear, logarithmic, and square- 
root law monochromatic limb darkening coefficients for 
a grid of Kurucz model atmospheres as published by 
IVan Hammd 11993). The Roche code, with the addi- 
tional input of the {u, v) coordinates, computes the vis- 
ibilities for a grid of wavelengths supplied by OYSTER, 
which are subsequently integrated over the NPOI band- 
passes. 

The fitting of the model parameters (except for the 
gravities) utilizes the Marquardt-Levenberg algorithm 
(jPress et al]ll992|) implemented in OYSTER, with the 
derivatives computed numerically. In addition to the 
visibility and phase measurements the reductions were 
also constrained to reproduce the observed V magnitude. 
This provided particularly strong constraints on the po- 
lar effective temperature. 

5. DISCUSSION 

5.1. Model Fitting 

The analysis proceded with few complications. In par- 
ticular, a close examination of the surface indicated 
no unusual morphology, and indeed the iterations con- 
verged to the same final solution independent of our 
starting guess, whether from larger or smaller values of 
the parameters. Our final solution, given in column 2 of 
Tabled is based on the triple phases and triple ampli- 
tudes only, the latter with an overall floating multiplier 
for each scan as described in i)3.6l The reduced for 
this solution is still a bit large, which we attribute to 
some residual non-gray problems of uncertain origin in 
the amplitudes, as shown in Figure^ The fit to the clo- 
sure phases, Figure |51 on the other hand is remarkable, 
showing no trends with wavelength or hour angle. 

The x^/^ for our best model is well in excess of unity. 
This is the result of both the remaining residuals in the 
triple amplitudes evident in Fig. ^ and the very small 
formal errors in both the amplitudes and phases owing 
to the high signal levels. If the large ha-d resulted 
from just the latter, we could have taken the usual expe- 
dient and scaled the errors by a common multiplier which 
would have normalized away the excess x^ • Estimates of 
the uncertainties in the resulting parameters would then 
have been obtained through the usual process of varying 
each parameter until x^ had increased by unity, thereby 
mapping out out formal errors and the correlation ma- 
trix. 



However, this gives very small estimated errors which 
we feel do not properly reflect the influence of the (small) 
remaining biases in the triple amplitudes. We have there- 
fore decided to take a conservative approach to our error 
estimates. We have run a separate reduction using the 
squared visibilities for the individual baselines in place 
of the triple amplitudes. This solution, which has a sub- 
stantially larger reduced x^ reflecting the larger residuals 
in the amplitudes of the individual baselines, is summa- 
rized in column 4 of Tabled Our adopted errors, shown 
in column 3, are the difference between these two solu- 
tions. 

Figure 01 shows how Altair appears projected on the 
sky. The intensity distribution at 500 nm, as would be 
seen for example by an interferometer, is color encoded: 
blue for high intensity, red for low intensity. Except for 
limb-darkening, this is also a temperature encoding. The 
range in intensities, a factor of 18, is about a factor of 
2.5 more than would be expected due to limb-darkening 
alone in a non-rotating star of this spectral type. 

Included in Table ^ is the integrated B — V a.s calcu- 
lated for the models. Although we force all models to re- 
produce the observed V magnitudes, as described above, 
this constraint does not automatically mean B — V will 
be reproduced. The fact that the color does agree with 
the measurements is a significant consistency check on 
the models, particularly given that the range of temper- 
atures across the surface could produce a wide range in 
B — V. A number of ancillary parameters derived from 
the adopted fit are given in Table El and are discussed 
further in the ^ 

5.2. Comparison with Previous Results 

As noted, Altair has been the subject of a number of in- 
terferometric measurements over the years, most notably 
by the PTI array (van Belle et al. 2001). The diameter 
determinations are not simply intercomparable, even if 
we consider only the projected major and minor axes of 
our model, since the PTI observations were fit to the 
physical dimensions of a Roche model, which included 
limb-darkening but not gravity darkening. Still, our de- 
termination of an equatorial radius of 1.988 i?0 seems in 
reasonable accord with their quoted 1.88 i?0. 

More problematic are the reported position angles of 
the rotational pole. First, w e note that in a preliminary 
report of this work ( Peterso n et al.ll2'004fl the pole is off 
by 180° due to a sign error. More complicated is the 
disagreement between the position angle from the PTI 
measurements, —25° ± 9 compared to our 123?2 ± 2?8. 
Communication with van Belle (2003, private commu- 
nication) indicated that the {u, v) coordintates of each 
baseline were inadvertantly exchanged in their analysis. 
Fitting a simple uniform ellipse model to the PTI data, 
corrected for the component swap, yields PA ~ 122?2, 
in excellant aggrement with our res ult. 

Recently Rcincrs & Roved ()2004fl have reported the 
determination of Altair's equatorial rotational velocity, 
Veq < 245 km s^^. They analyzed the star's rotational 
broadening profile to determine the first two zeros of its 
Fourier transform, the ratio of which has been shown 
by Rcincrs (200^ to depend on the equatorial velocity 
rather than the usual vsini. This is a new approach to 
measuring total velocities in stars, and it is difficult to 
know how much weight it should be given. One notable 
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TABLE 1 

Roche Model Fits for Altair and ( Aql 







Altair 






C Aql 




Parameter 


Vl23 & 4>cl'' 


Errors V^'' 


&^<t>ci'' 13 = 


0.09 


V," & (/-cl"^ 


Errors 


Lu = n/Qc 


0.90 


±0.02 


0.88 


0.978 


0.990 


±0.005 


6p (mas) 


2.96 


0.04 


3.00 


3.04 


0.815 


0.005 


Tp (K) 


8740 


140 


8600 


7980 


11750 




I (deg) 


63.9 


1.7 


62.2 


65.6 


90 


+0, -5 


PA (deg) 


123.2 


2.8 


120.4 


97.4 


45 


±5 


X^/DOF 


3.8 




9.5 


13.4 






V (obs: 0.77) 


0.765 




0.765 


0.76 






B-V (obs: 0.22) 


0.215 




0.22 


0.26 






Veq sini (kms~-^) 


245 




231 


295 







^Model fit to the triple amplitude and closure phase data. 

''Models fit to the closure phase data and the squared visibilities of the three baselines. 



aspect of that analysis was the adoption of an exponent 
for the gravity darkening law (0 ~ 0.09) that was about 
1/3 that of the lvon Zeipell l)1924D value. We discuss this 
aspect of calculating rotationally distorted stars next. 

5.3. Gravity Darkening 

While we have several lines of reasoning, described in 
ti4.1l that lead us to believe that solid body rotation is 



valid for these models, the situation is not so clear with 
regard to the exponent on the gravity darkening law, 
Teff oc g^ff The classical work of Ivon Zeipell l)1924f) 
would seem to apply to an A star, even one somewhat 
evolved, yielding /3 — 0.25. However, convection does 
occur in late A stars, and particularly in the photosphere 
where it competes with radiation in carrying the flux. It 
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is then important to note that ILucvl l|1967j) has shown 
that for small distortions the appropriate coefficient in 
fully convective envelopes is closer to /3 ~ 0.08. 

This leaves matters in a somewhat uncertain state. On 
the one hand, there are stars in the transition region 
between having fully convective and fully radiative en- 
velopes, Altair arguably one of them, and there is no 
obvious guidance in choosing an appropriate value for 
this parameter. On the other hand, even where the en- 
velopes are unambigiously in either one of those states 
or the other, the classical results apply under rather dif- 
ferent circu mstances, which w e nex t discuss. 

Both the Ivon Zeipell 1)19241) and ILucvl l|1967|) results 
treat rotation as a perturbation. However, in the ra- 
diative case where uniform rotation is adequate (and is- 
sues like mass loss, etc. can be ignored), the quantity 
treated as a perturbation is the size of the quadrapole 
moment of the gravitational field. Since stars are cen- 
trally condensed, even for velocities approaching critical 
the distortions in the core are modest, and one can ex- 
pect that the analysis given bv lvon Zeipell l|1924r ) will be 
reasonab ly accurate. Thi s has been found to be true in 
practice llSackmannlll"97(il) . 

For the convective case, the gravity darkening expo- 
nent is obtained by analyzing the adiabats fou nd in the 
envelopes of representative stars. ILucvl l)1967l) quite ex- 



plicitly points out that the derivation is valid only for 
small changes in the effective gravity. Of course, the 
effective gravity changes by orders of magnitude as rota- 
tion approaches critical, and it is not clear whether the 
exponent derived by Lucy can be used to describe grav- 
ity darkening for anything but the most modest rotation. 
Again, this is in contrast to the modest contribution of an 
induced gravitational quadrapole, even for stars rotating 
at breakup. 

Even so, in a recent series of papers, IClaretl ll200l and 
references therein) has attempted to deal with the is- 
sue of a smooth interpolation between these two extreme 
cases. He has noted that as stars evolve off the main 
sequence and toward the red giant branch, their interior 
structures trace out approximately straight line loci in 
a (log Te//, log 5) diagram. On the main sequence for 
massive (mostly radiative) stars the slope of this line is 
about 0.25, and for intermediate mass stars (~ lAf©, 
mostly convective) the slope is about 0.06, the two val- 
ues being remarkably close to the radiative and convec- 
tive e xponents cited above. Working with the interiors 
codes, IClaretl l)2004D is able to evaluate this exponent at 
each point in the evolutionary paths of models covering 
40 > Mq > 0.08, offering the results as appropriate ex- 
ponents to use in rotating stars and stars in close binary 
systems over their entire evolutionary lifetimes. 



Altair i=63.9 



Resolving Altair 9 
25-May-01 AE-AW Baseline 




RA (mas) 

Fig. 3. — A false-color rendering of Altair's visible surface. In- 
tensity at 500 nm increases from red to blue. Except for the effects 
of limb-darkening, this is also a map of temperature, which varies 
from 8740 K at the pole to 6890 K at the equator. 



This is a constructive suggestion for the thorny prob- 
lem of choosing an appropriate gravity darkening law. 
However, we are not fully convinced of the leap of going 
from deriving a quantity based on evolutionary changes 
to using it to describe the effects of rotational distor- 
tion. One might make the case for small amounts of ro- 
tation or if it could be shown that rotational distortions 
and evolutionary effects were close to being homologous 
transformations from one to the other. But rotational 
distortions are not homologous to evolutionary changes 
and it is not at all clear how well these "interpolations" 
work. 

The observations reported here bear somewhat on this 
problem. We have tried converging o ur Roche models 
using the value P = 0.09 adopted by iReiners fc Roved 
(|200l . The results are shown in column 5 of Table 
The data set was used so the comparison is with 
column 4. To achieve the degree of asymmetry found 
in the triple phase data with this low exponent value, 
the rotation parameter is forced to near critical rotation, 
CO ~ 0.978. In turn the predicted projected rotational 
velocity, vsini ^ 295kms~^ conflicts with the observed 
value, and the predicted color is significantly redder than 
observed. Further, the reduced is significantly worse 
for this fit. 

We feel it is premature to use these observations to 
derive a "best" value of the gravity darkening parameter 
until the remaining visibility residuals are better under- 
stood. Since it is the phase measurements, which we 
trust, which are sensitive to asymmetries in brightness 
across the disk and not the visibility amplitudes, it does 
appear that the von ZciDcl (1924) value for that param- 
eter is superior to the value adopted bv IReiners fc Roverl 
ll2004) . 

Moreover, as we will describe in ti6.2l our results indi- 
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Fig. 4. — One of the observed V^'s for the AE-AW baseline. 
The effect of the strong H/3 feature in the A 468.3 channel is clearly 
evident and well matched by the model. 



cate that Altair has hardlv evolved from the zero age 
main sequence. The iClareti lll998l l2004|) tables give 
values for /3 in essential agreement with the Fvon Zcioej 
l)1924ft result during this phase of evolution. Thus, Al- 
tair does not provide a test of those tables, but does 
highl ight a difficulty i n app lying the technique proposed 
by IReiners fc Roved 1)20041) for finding total velocities, 
namely the difficulty of obtaining a priori reliable esti- 
mates for f3. 

5.4. Hf3 

In Figure 01 we show the blue squared visibilities of 
the AE-AW baseline for the HA = 0'.'83 observation. 
The notable feature at 486 nm is H/3; the agreement with 
the calculations shown in this scan, and the others not 
shown, is striking. This feature is nearly centered in the 
486.3 nm channel. In contrast. Ha, which has a smaller 
equivalent width and is split between channels at A 665.4 
and A 648.7, is much less noticable. The main reason for 
the feature being reflected in the amplitudes is the reduc- 
tion in the limb-darkening coefficient, the star appears to 
be more like a uniform disk at this wavelength, and thus 
the visibility is reduced. This close agreement is a nice 
confirmation of the details of the model fits. 

6. CONCLUSIONS 
6.1. Imaging Altair 

The primary result of these observations is, we be- 
lieve, the first detection of asymmetric surface intensi- 
ties on the surface of a star induced by rotation. Al- 
though we have imposed a model on the data and fit the 
m odel parameters, the simple conclusion, first reported 
bv lOhishi et al.l l|2004D . is that the surface of Altair dis- 
plays an extremely asymmetric intensity distribution and 
that the asymmetry is consistent with that expected from 
the known high rotation and with the previously reported 
oblateness ijvan Belle et al.l200lD . We have, in effect, im- 
aged the surface of an A star. 

In Table El we summarize various physical parameters 
for the adopted model — column 2 of Table H Most 
quantities should be self-explanatory. Subscript "B" 
refers to the model if it were rotating at breakup. The 
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angular diameter 0min is for the projected minor axis 
while 9 Max is for the projected major axis, i.e., the an- 
gular diameter of the equator. 

6.2. The Status of Altair 

Knowing now the rotational state of Altair, we can 
better answer fundamental questions such as its evolu- 
tionary status. Over the years Altair has been classi- 
fied as A7 IV-V (e.g., Johnson & Morgan 1953), the lu- 
minosity class usually indicating an object slightly past 
the end of its main sequence evolution while the spectral 
type is that of a star having an effective temperature 
in the vicinity of 780 0K (e.g., lErspamer fc Northll200l 
Ivan Belle et al.ir2001|) . In the context of analyzing Al- 
tair's pulsations, it is important to know in detail its 
evolutionary state, mainly the extent of its core, and to 
remove any biases that might be introduced by rotation. 

Fortunately, in the context of rigid rotat ion, th is is not 
so difficult. Early results fe.g.. ISackmannl [l97Cl() showed 
that two quantities were relatively insensitive to the ef- 
fects of rotation: polar radius and total luminosity. The 
Roche model fits give polar radius directly while it is a 
strai ghtforward m atter to calculate the total luminosity 
fe.g.. lCollin3 ll963h. 

These two quantities are not perfectly conserved. In 
the range of interest, 3 Mq > M > 1.4 M©, both quan- 
tities decrease with increasing rotation, approximately 
in proportion to w^, reaching a maximum co rrection of 
abou t 6% in luminosity and 1.5% in radius ijSackmannI 
For stars in the neighborhood of 1.8 Afp, rotating 
with ~ 0.8 we find from the Sackmann ( 1970) cal- 
culations that the non-rotating star would be 4% more 
luminous and 1% larger than our deduced polar radius. 

To estimate the parameters of the appropriate non- 
rotating star we have used the evolutionary tables by the 
Geneva group (Schaller ct al. 1992). These models were 
calculated with modest convective overshoot (0.2iJp). 
We have used the grid for a composition oi X — 0.68, 
Y = 0.30, Z = 0.02. 

The quantities given in Table El as "Non-rotating" are 
those estimated from the Schaller ct al. (1992) models 
and are quite striking. The last entry in Table [21 is the 
mass fraction of hydrogen remaining in the core. This is 
to be compared to a starting value of = 0.68. Altair 
is almost on the ZAMS. 

6.3. Chromospheric indicators 

As mentioned in Altair is one of two A7 objects, 
the other being a Cep, in which certain ultraviolet emis- 
sion lines, taken as indicators of chromospheric temper- 
ature inversions, are seen. No objects of earlier spectral 
type are known to show these features, and it is usually 
argued that these therefore represent the hottest photo- 
spheres where convection is still capable of creating such 
temperature profiles. The model for Altair adopted here 
calls that conclusion into question. As shown in Table 
12 Altair has a broad swath of 6900 K gas at its equa- 
tor, which is the likely source of the str ong convection. 
We a lso note the recent announcement l)van Belle et al.l 
12005ft that substantial oblateness and gravity darkening 
have been found in a Cep, suggesting significant amounts 
of cool, convective gas in that object as well. 

6.4. A S Scuti star 



TABLE 2 

Altair Physical Parameters'' 



Quantity Unit Value Error 



Rotating Model Parameters 





kms ^ 


273 


13 




km 


374 


3 


n 


day~i 


2.71 


nil 




day~i 


3.01 


0.06 


Tp 


K 


8740 


140 


Teg 


K 


6890 


60 


Rp 


Rq 


1.636 


0.022 


Reg 


Re 


1.988 


0.009 


^min 


mas 


3.056 


0.047 




mas 


3.598 


0.017 


logL 


Lq 


1.027 


0.011 




cgs 


4.266 


0.012 


loggeq 


cgs 


3.851 


0.035 


Non-rotating Parameters 


C 


M 


Mq 


1.791 


0.018 


R 


Rq 


1.652 


0.022 


log 9 


Cgs 


4.256 


0.002 




K 


8200 


98 






0.607 


0.019 



^Uncertainties due to the parallax 
have not been included in the errors. 

''Rotating at breakup but with the 
same mass and polar radius. 

'^The parameters of a non- 
rotating star fr om t he Geneva grid 
I Schaller ct al.' 119921) which would 
reproduce the (corrected) luminosity 
and polar radius - see text. 



As we also indicated in ^ rather than just being 
another rotating A star, Altair may prove to be a 
valuable laboratory for examining the internal rotation 
state of a star w i th a predominantly radiative enve- 
lope. iBuzasi et al.l l)2005tl have announced the discovery 
of S Scuti pulsations in Altair and have identified several 
of the periods with frequencies mostly in the range 15 - 
18 day"-'^. Interestingly, two of the frequencies reported 
by Buzasi et al. (2005), 3.526 day^^ and 2.57day~^ were 
significantly lower, the latter being quite close to the ro- 
tational frequency we have derived, e.g., Table (21 

This was immediately recognized as providing a poten- 
tial probe of the interior structure and particularly the 
rotation law, and an atte mpt has been mad e to iden- 
tify and model the modes ijSuarez et al. 2005) . Unfortu- 
nately, not knowing the rotational state of the star and 
making the assu mption that equator- on was the most 
likely orientation, ISuarez et alJ lj2005f) adopted a total 
rotation significantly lower than now seems likely. Other 
effects of this choice included identifying the evolution- 
ary state as being substantially more advanced than we 
believe is the case. As is clear from their results, rota- 
tional velocities above the 180-240 km s~^ range they in- 
vestigated lead to large changes in the oscillation modes, 
making mode identification difficult. It will not be an 
easy task to tap the information being provided by Al- 
tair. 
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The referee, Gerard van Belle, asked a number of pen- 
etrating questions that significantly improved this work. 
Notably, he wondered about the effects of rotation on our 
calibrators, leading to the discovery that we could mea- 
sure those effects in C, Aql. The NPOI facility is a collab- 
oration between the Naval Research Laboratory and the 



US Naval Observatory in association with Lowell Obser- 
vatory, and was funded by the Office of Naval Research 
and the Oceanographer of the Navy. This research has 
made use of the SIMBAD literature database, operated 
at CDS, Strasbourg, France, and of NASA's Astrophysics 
Data System. 
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